Strain Stiffening Induced by Molecular Motors in Active Crosslinked 
Biopolymer Networks 
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We have studied the elastic response of actin networks with 
both compliant and rigid crosslinks by modeling molecu- 
lar motors as force dipoles. Our finite element simulations 
show that for compliant crosslinkers such as filamin A, the 
network can be stiffened by two orders of magnitude while 
stiffening achieved with incompliant linkers such as scruin 
is significantly smaller, typically a factor of two, in excel- 
lent agreement with recent experiments. We show that the 
differences arise from the fact that the motors are able to 
stretch the compliant crosslinks to the fullest possible ex- 
tent, which in turn causes to the deformation of the fila- 
ments. With increasing applied strain, the filaments fur- 
ther deform leading to a stiffened elastic response. When 
the crosslinks are incompliant, the contractile forces due to 
motors do not alter the network morphology in a significant 
manner and hence only small stiffening is observed. 

The mechanical properties of plant and animal cells are gov- 
erned by the cytoskeleton, a flexible and dynamic network of 
biopolymer fibers combined with a group of associated regula- 
tory and crosslinking proteins 1 2 . One of the key aspects of the 
mechanical behaviour of these networks is their highly non- 
linear elastic response to applied stresses^, in particular their 
ability to strain stiffen by orders of magnitude when subject 
to large stresses. Cells also employ molecular motors to con- 
vert chemical energy into mechanical work 1 . Motors generate 
internal stress in the networks even in the absence of external 
loadingEE] In this manner, cells can regulate their mechanical 
properties by using both active and passive components. 

While the mechanical behavior of semiflexible polymer net- 
works with compliant and rigid crosslinks has been studied in 
detail both experimentally and theoretically^"^!, the interplay 
between active mechanisms of stress generation through motor 
activity and passive strain hardening properties of crosslinks 
has only been considered very recently 47 . In this regard, re- 
constituted actin networks can be particularly useful, since the 
density of crosslinks and motors can be varied in a desired man- 
ner to gain insights into the mechanisms of strain hardening and 
nonlinear elastic response. Indeed, recent experiments on net- 
works that consist of actin filaments crosslinked by filamin A 
(FLNa) and bipolar filaments of muscle myosin II show that in 
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the absence of any applied loads 7 , the motors stiffen the net- 
work by about two orders of magnitude. The degree of stiff- 
ening was found to increase with increasing density of myosin 
motors. Another key observation from this study relates to the 
magnitude of stiffening caused by compliant and incompliant 
crosslinks. While FLNa is a large, highly flexible dimer that 
promotes orthogonal F-actin crosslinking, scruin is an incom- 
pliant crosslink. Interestingly, it was found that in distinct con- 
trast to FLNa, scruin does not promote active stiffening of F- 
actin networks upon addition of myosin. These results clearly 
show that actomyosin contractility when combined with appro- 
priate crosslinks can allow the cell to operate in a nonlinear 
regime to actively control its mechanical response. 

Why do compliant crosslinks in active networks lead to large 
strain stiffening while no significant increase in stiffness is 
observed in the case of incompliant crosslinks? To answer 
this question and to quantitatively study the interplay of inter- 
nal strains generated by molecular motors and external loads, 
we study the mechanical response of active networks using fi- 
nite element simulations where the motors are treated as force 
dipoles (Fig. [TJ. A number of approaches including mean- 
field models, effective medium theory and numerical simula- 
tions have been used to study the elastic response of passive 
networks with both compliant and incompliant crosslinks^. 
However the role of crosslinks on the mechanical response of 
networks with molecular motors has not been considered the- 
oretically. Our work shows that the nature and density of the 
crosslinks play a key role in determining the strain stiffening re- 
sponse in active biopolymer networks. FLNa is a crosslink that 
is compliant at small pulling forces, but is stiff beyond a critical 
value of stretch 1 1 . We find that even in the absence of applied 
loads, motors lead to almost completely stretch out the com- 
pliant crosslinks taking them into the stiffened regime, which 
in turn also leads to bending of the filaments. As the crosslinks 
are fully stretched, the applied load is accommodated by the de- 
formation of filaments. Since it is more difficult to deform the 
filaments compared to stretching of compliant crosslinks, the 
network stiffens by as much as two orders of magnitude com- 
pared to the case when motors are absent. On the other hand, 
when the crosslinks are incompliant, the contractile forces due 
to motors do not alter the network morphology in a significant 
manner leading to much lesser stiffening (typically by a factor 
of two) in agreement with experiments^. 

The biopolymer networks considered in our study are con- 
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Figure 1 Schematic of actin filaments (black lines) with myosin mo- 
tors (blue) and compliant crosslinks (green lines). Myosin motor ex- 
erts equal and opposite forces on the filaments on which it is attached, 
which results in a force dipole (pair of red arrows). These forces lead 
to the extension of the filamin crosslinks. 



structed as follows. Straight filaments of length L = 10 fjm 
(comparable to the actin filament lengths of ^10-15 /jm in 
Ref.^ and random orientations are assembled in a square box 
of width W = 40 /jm. When two filaments cross each other, 
they are connected by a nonlinear spring (to be described be- 
low) or are linked rigidly. Filaments that cross the top and 
bottom boundary are cut and the dangling ends are removed, 
while periodic boundary conditions are imposed along the lat- 
eral boundaries. The system is loaded by restricting the hor- 
izontal displacements of the nodes of the filaments at the top 
and bottom of the box to yW and zero, respectively, where y is 
the applied shear strain. We assume that the polymer filaments 
that make up our network are semiflexible, so that the persis- 
tence length of the individual chains is much longer than the 
average distance between the two crosslink sites, and compara- 
ble to the contour length of individual polymer chains. There- 
fore we ignore thermal energy arising from fluctuation of the 
filaments and consider only the extensional and bending ener- 
gies of the filaments. Using K and /j to denote bending modu- 
lus and stretching modulus of F-actin, respectively, we choose 
fi/L = 1.6 MPa and ^k/(/jL 2 ) = 2.3 x 10~ 4 (both represen- 
tative of F-actin networksU^. The normalized line density of 
our network (p = LtL/W 2 , where Lj is the total length of fila- 
ments in the cell) is p = 12.5, which is well above the rigidity 
percolation threshold of p = 5.7^= 

The contour length of FLNa is /o = 150 nm and its per- 
sistence length is l p = 20 nmPl. For stretches less than the 
contour length, this crosslink behaves like a linear spring with 
spring constant k c i = 13 , where kgT = 4.11 pN ■ nm. 
When stretched beyond its contour length Iq, the stiffness of 



the crosslink increases very rapidlyEE] Following the work of 
Brodersz et al. 6 , we therefore model the crosslinks as piece- 
wise linear springs such that the force F = —kAl, where k = k c i 
for length Al < Iq, while k S> k c i when Al > Iq, where Al is 
the extension of the crosslink. Scruin crosslinks on the other 
hand are inextensible 3 . Very recent measurement shows that 
force applied by a single skeletal myosin head on an F-actin is 
in the range lpN - 5pN 14 . Myosin is typically assembled into 
thick filaments with several heads; the filaments typically have 
lengths of ~ IfmP^. To model the contractile force gener- 
ated by the motors we model them as force dipoles as shown 
schematically in Fig. [T] In our simulations, the motors are as- 
sembled by picking a point at random on a filament and then 
picking another point located on a neighboring filament such 
that the distance of the two points lies in the range of 0.5-2 /jm. 
A typical force of lOpN is applied along the line connecting the 
two points to simulate the pulling effect of the motors as shown 
in Fig. [1] 
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Figure 2 Stress (f) and differential shear stiffness (K') versus shear 
strain (y) for active F-actin network with rigid or compliant crosslinks. 
The myosin motors only stiffen the F-actin network with rigid 
crosslinks by a factor of 1.4 at small strains: network without myosin 
motor (black dashed line) and N m /N act i„ = 6.4 (green circles). How- 
ever the myosin motors stiffen the flexibly crosslinked F-actin network 
up to a factor of 70: network without myosin motors (black solid line), 
with N m /Nactin = 0.9 (red solid line), with N m /N ac iin = 3.1 (blue solid 
line), and with N m /Nactin = 6-4 (green solid line). N m and N ac tin corre- 
spond to the total number of motors and filaments, respectively. Each 
motor exerts a force of lOpN on the actin filament. 

The elastic fields in the sheared filament network are com- 
puted using the finite element method, discretizing each fila- 
ment with 100 equal-sized Timoshenko beam elements^! All 
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simulations are carried out in a finite deformation setting; i.e., 
the effect of geometry changes on force balance and rigid body 
rotations are explicitly taken into account. The macroscopic 
shear stress X for our model is the total horizontal reaction 
force at the top of the simulation box divided by W. The di- 
mensionless stress X in Fig. [2] is defined as T = xL/fi, and the 
modulus K' = dx/dy. To gain insight into different deforma- 
tion mechanisms, we also compute the total energy associated 
with stretching of the crosslinks in the cell, E c , the total bend- 
ing energy of all the filaments, Eh, and the sum of the bend- 
ing and stretching energies or the total strain energy of all the 
filaments, Et. We first consider the case where there are no 
motors present and the crosslinks are either rigid (black dashed 
lines in Figs.[2]and[3]l or compliant (black solid lines in Figs. [2] 
and^. As in earlier work 6 , the modulus in the latter case shows 
a sharp increase when the crosslinks are fully stretched and en- 
ter the hardened regime. It can be seen from Fig. [3]that below 
this threshold the total energy of the network is dominated by 
the energy of stretching the crosslinks, but above this threshold 
the filaments first bend and then stretch. In the case of rigid 
crosslinks, for small strains, the deformation of the network is 
dominated by bending of the filaments followed by stretching 
and orientation of the filaments along the direction of shearing. 
Note that at very large applied strains, both the networks show 
identical response due to stretching of the filaments. 

Next we consider the response of the network with compli- 
ant crosslinks in the presence of motors. At small strains, Fig. [2] 
shows that the networks continuously stiffen with increasing 
density of motors. Even with one motor per filament, the stiff- 
ness of the network increases by over one order of magnitude. 
Networks with about six motors per filament are stiffer than the 
networks without motors by about two orders of magnitude. To 
understand the reason for this marked increase in stiffness, we 
consider a) the ratios of energies associated with different de- 
formation modes (Fig. [3}, b) the distributions of the lengths of 
the crosslinks (Fig. HI and c) the shapes of the deformed fila- 
ments (Pig. |5j. Fig.[3]shows that even at small applied strains, 
the ratio of the total energy of the crosslinks to the total strain 
energy of all the filaments decreases significantly with increas- 
ing motor density. Indeed, the histogram of crosslink exten- 
sional lengths in Fig. HJ confirms that the motors are able to 
induce stretching of almost all the crosslinks to their contour 
length, k), beyond which it is difficult to stretch them. This 
causes bending of the filaments, which explains the decrease in 
the ratio of the total energy of stretching of crosslinks to the to- 
tal energy of the filaments, E c /Ef, with increase in the density 
of motors. Upon application of a small external load, the defor- 
mation is primarily borne by the bending of the filaments lead- 
ing to stiffer response compared to the case where motors are 
absent. The dominance of the bending modes is clearly seen in 
Fig [5] for small applied strain of 0.05, the filaments in the net- 
work without motors are straight whereas significant bending 



(a) 


10 1 




10° 




10" 1 




10 -2 


(b) 


10 -3 


10° 




10 -1 




ltr 2 




1(T 3 




ltr 4 

















a 1 

T 
\ 

\ 



10 



10 



10" 



Figure 3 (a) The ratio of total stretching energy (E c ) of compliant 
crosslinks to total energy (Ef) of filaments in the system, and (b) the 
ratio of total bending energy (£7,) of the filaments in the system to to- 
tal strain energy of all the filaments (EA as a function of shear strain 
(y). At small strains, the rigidly crosslinked F-actin network is domi- 
nated by bending of F-actin, regardless of whether myosin motors are 
present (green circles) or absent (black dashed line). However, myosin 
motors drive the deformation from stretching of the crosslinks to bend- 
ing of F-actin as evidenced by the decrease in the ratio E c /Ej from the 
case with no motors (black solid line), to the cases with increasing con- 
tractile forces (not plotted here) or motor densities (N m /N ac ,i„ = 0.9 
(red solid line), N m /N a ain = 3.1 (blue solid line), and N m /N act i n = 6.4 
(green solid line)). All flexibly crosslinked F-actin networks with or 
without myosin motors show nearly identical behavior of the ratio 
EbjEf (green solid line in (b)). 



of the filaments can be seen in the network with motors. Also 
note from Fig.HJthat most of the crosslinks in the former case 
have not been stretched to the maximum extent. As the applied 
strain is large (~ 0.3), in all cases there is a transition in the de- 
formation modes of the filaments from bending to stretching, at 
which point the response of all the networks is identical. 

In contrast to the filaments with compliant crosslinks, the 
stress-strain curves and the incremental moduli of networks 
with rigid crosslinks are not significantly altered by the pres- 
ence of motors. Indeed, we find only an increase close to a 
factor of two in the modulus of the network at small strains. 
We can understand this by noting that the filaments in rigidly 
crosslinked networks deform primarily through bending. The 
motors do not qualitatively change this picture. The effect of 
the motors on bending deformation is small unlike the large 
change in the lengths of the crosslinks that they induce in the 
case of compliant networks. Our calculations therefore pro- 
vide an explanation for the key role played by the nature of the 
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Figure 4 Distribution of the extension of the crosslinks Al relative 
to their contour length Iq at shear strains y = (blue) and y = 0.05 
(green). N c i is the number of crosslinks normalized by the total num- 
ber of crosslinks in the active network. The symbols A-C corre- 
sponded to the black squares given in Fig. [2] In B and C myosin mo- 
tors drive most compliant crosslinks up to contour lengths, /q, whereas 
when no motors are present (A), the crosslinks have not been stretched 
to their fullest extent. Note that more crosslinks are stretched to their 
contour lengths in C than in B owing to the larger motor density in the 
former case. 



crosslinks on the mechanical response of active networks. 
a c D 




Figure 5 The deformation of F-actin network corresponding to the 
points A, C, and D in Fig. [2] Deformation is dominated by stretching 
of compliant crosslinks (A), bending of flexibly crosslinked F-actin 
(C), or bending of rigidly crosslinked F-actin (D). 

It is well known that in uncrosslinked networks the myosin 
filaments effectively fluidize actin networks by actively sliding 
antiparallel actin filaments past one another, which can lead 
to large scale reorganization of the network. However in a 
crosslinked network, upon addition of myosin, there was no 
noticeable change in network structure 7 . When myosin and 
FLNa were both present, the network still remained homoge- 
neous and unbundled suggesting that relative sliding between 
the motors and filaments is relatively small, which justifies the 
treatment of the motors as force dipoles. It is also possible that 
myosin-driven tension may release FLNa crosslinking. This is- 
sue, the effect of potential sliding and perhaps the stretching 
of myosin motors themselves can be considered in future work 
using the model we have developed here. 

In summary, by using material parameters typical for actin 
networks with compliant crosslinks, we have shown that 
myosin II motors generate internal stresses by stretching the 
crosslinks which in turn pull on the actin filaments. Once the 



crosslinks are fully stretched to their contour lengths, the dif- 
ferential stiffness of the network can increase by two orders 
of magnitude in excellent agreement with recent experiments^. 
In addition, our simulations show that motors do not lead to 
any significant stiffening in the response of rigidly crosslinked 
networks, also in accord with experiments 7 . These observa- 
tions underscore the importance of the nature of crosslinks 
on determining the strain hardening behavior of active net- 
works and provide guidelines for tuning mechanical response 
of biomimetic systems. 
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